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■ Abstract 
» I 

O , Binary black hole spacetimes with a helical Killing vector, which are discussed as 



an approximation for the early stage of a binary system, are studied in a projection 



o 

lO ' formalism. In this setting the four dimensional Einstein equations are equivalent to a 

. three dimensional gravitational theory with a SL{2, C)/SO{l, 1) sigma model as the 

I material source. The sigma model is determined by a complex Ernst equation. 2 + 1 

' decompositions of the 3-metric are used to establish the field equations on the orbit 

■ space of the Killing vector. The two Killing horizons of spherical topology which 

! characterize the black holes, the cylinder of light where the Killing vector changes 

I from timelike to spacelike, and infinity are singular points of the equations. The 

horizon and the light cylinder are shown to be regular singularities, i.e. the metric 
, functions can be expanded in a formal power series in the vicinity. The behavior 

Q I of the metric at spatial infinity is studied in terms of formal series solutions to the 

O^' linearized Einstein equations. It is shown that the spacetime is not asymptotically 

Ch^, flat in the strong sense to have a smooth null infinity under the assumption that 

the metric tends asymptotically to the Minkowski metric. In this case the metric 
functions have an oscillatory behavior in the radial coordinate in a non-axisymmetric 
^ . setting, the asymptotic multipoles are not defined. The asymptotic behavior of the 

^ I Weyl tensor near infinity shows that there is no smooth null infinity. 



1 Introduction 



Binary black hole systems in the last stage before coalescence are the most promising 
sources of gravitational radiation to be detected with the first generation of gravitational 
wave detectors. From a theoretical point of view this is a difficult relativistic problem 
which can possibly only be solved numerically since there are no symmetries. The ad- 
vantage of black hole systems is that no matter is involved and that only the vacuum 
equations have to be studied. It is generally expected that binary systems will have an 
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early stage of quasi-circular motion for large separation of the binaries. Radiation damp- 
ing will lead to almost circular orbits with the radius decreasing in time due to the emitted 
radiation. For smaller distances of the binaries this motion is expected to be followed by 
a rapid inspiral. The end result will be a single black hole which will settle to a stationary 
hole in a ring down phase which can be described by black hole perturbation theory. For a 
given binary system an important characteristic quantity is the innermost stable circular 
orbit (ISCO), the last almost circular orbit before the final inspiral. In general there is 
no unique definition of the ISCO, but it should be possible to define a characteristic scale 
where the quasi-stationary early phase of the binary system comes to an end. 

In the case of the binary motion of two oppositely charged particles, Schonberg pP and 
Schild |2] considered an approximation to the quasi-stationary phase of the system. The 
quasi-circular motion is approximated by a sequence of exactly circular orbits which are 
obtained by exactly compensating the outgoing radiation by ingoing radiation. The bind- 
ing energy of the system as defined in decreases with the distance of the charges up to 
some minimal value which can be taken as the definition of the ISCO: for smaller values 
of the distance, more and more incoming radiation is needed to stabilize the circular mo- 
tion. The approximation thus predicts its own breakdown and allows for an unambiguous 
definition of the ISCO. The quasi-stationary approximation corresponds to a so-called 
helical Killing vector ^ of the system which is in standard Minkowski coordinates given 
by ^ = 9i + ^d^. The main features of such a vector can already be inferred from this 
case: the vector becomes null at the so-called light cylinder given by p = 1/Q where the 
observer rotating with the angular velocity Q rotates with the velocity of light. In the 
interior of the light cylinder, the Killing vector is timelike, in the exterior spacelike. In a 
general spacetime with a helical Killing vector, the light cylinder will be deformed, but 
will still have cylindrical topology. 

Detweiler [SI IH E] suggested to use this concept to describe the quasi-circular regime of 
binary black holes which corresponds to studying spacetimes with a helical Killing vector. 
Since Einstein's theory is a nonlinear theory, the incoming radiation will lead to a space- 
time which is not asymptotically flat in a strong sense (mass and angular momentum can- 
not be defined in the usual way asymptotically). It was shown in |Bj that spacetimes with 
a helical Killing vector cannot have a smooth null infinity if there is no additional station- 
ary Killing vector close to X, see also Though the Arnowitt-Deser-Misner(ADM)-mass 
cannot be defined, Friedman et al. [S] could show that a thermodynamical treatment as in 
the single black hole case is possible and that there exists a first law. With the help of the 
first law, the ISCO can be defined for asymptotically flat spacetimes (which are possible 
for instance in the case of the first post-Newtonian order) as in the Maxwell the 
minimum of the binding energy which also marks the onset of dynamical instability. It 
is not clear how this result can be generalized to non-asymptotically fiat spacetimes in 
full general relativity. An important conceptional advantage of spacetimes with a helical 
Killing vector is the fact shown in |H| that spatially compact Killing horizons are event 
horizons. This allows for a local characterization of the event horizons in these models. 
It is thus not necessary to use local concepts as an apparent horizon or the concepts 
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developed in 0. 

In the study of binary black hole system, mainly approximative and numerical methods 
have been used so far. Post-Newtonian calculations have been carried out up to the 
third post-Newtonian order including resummation techniques (see [TIH ITU IT^ IT^ t^^- 
Within this approximation the determination of the ISCO appears to be self-consistent 
since corrections to its value due to higher order terms can be shown to be negligible. 
The post-Newtonian metric, however, cannot be used close to the horizons since black 
holes are the strongest relativistic objects known. Numerical calculations so far have 
been mainly performed within the Isenberg- Wilson-Matthews (IWM) theory ^3 , an 
alternative theory of gravitation without radiation. It follows from the Einstein equa- 
tions in a standard 3 + 1-splitting for a conformally flat spatial metric on the t = const 
hypersurfaces. Only the trace of the 6 time evolution equations is considered. By con- 
struction the theory coincides with the Einstein theory for spacetimes with conformally 
fiat spatial slices and thus reproduces exactly the first post-Newtonian approximation. It 
has to be noted that the Kerr solution does not allow conformally fiat spatial hypersur- 
faces, see [T7j. In 11111111120] initial values were constructed from the constraint equations 
via so called Bowen-York initial data by solving the Lichnerowicz equation numerically. 
This resulted in a significant discrepancy with post-Newtonian results for the ISCO. In 
PH 122] , complete IWM binary black hole spacetimes with a helical Killing vector were 
constructed numerically. The results were in good agreement with post-Newtonian re- 
sults, but suffered from an inconsistency in the model in the form of non-regular horizons. 
This non-regularity appears to be unavoidable for IWM spacetimes. To answer the ques- 
tion whether the Einstein equations allow for binary black hole spacetimes with smooth 
disconnected horizons in the presence of a helical Killing vector, it therefore appears nec- 
essary to study the fully relativistic situation. In !23j the constraint equations are solved 
in the presence of an approximate helical Killing vector on a background which is the 
superposition of two Kerr-Schild metrics. 

In general relativity, a helical Killing vector will lead to Einstein equations which are 
of a mixed type, i.e. elliptic in the interior of the light cylinder and hyperbolic in the 
exterior. Such equations are studied in aerodynamics ^24j in the case of transonic flows. 
For a review on mixed problems arising in gravitation, see [2n]- It was shown by Torre 
[2n] that the resulting equations can be written in the form of a first order differential 
system which belongs to the symmetric positive systems of Friedrichs |2Z| and Lax and 
Phillips j2H] • Classes of boundary conditions compatible with these equations are given 
in [2n]- Numerical studies of two-dimensional equations of this type were performed in 
|29| EO] • Three-dimensional toy models for the helically reduced Einstein equations were 
considered numerically in |31j . 

The purpose of this paper is to study the Einstein equations in the presence of a helical 
Killing vector for a vacuum spacetime with two disconnected Killing horizons of spherical 
topology. The idea is to derive a set of equations which is well suited for a numerical 
treatment by taking full advantage of the Killing vector. In the vicinity of the critical 
points of the equations as the black hole horizons, the metric functions will be given in 
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terms of formal power series which should be useful for the numerical implementation. 
The approach is similar as in the study of cosmological singularities, see e.g. [321 
references therein. These methods are purely local and do not imply an existence proof 
for a spacetime with a helical Killing vector and two regular horizons. The use of so-called 
Fuchsian methods (see [SS]) in the context of cosmological singularities would not change 
the situation since even an existence proof for one regular horizon would not imply the 
existence of the second regular horizon one is interested in here. Therefore we will not 
discuss the question of existence and convergence radii of the formal series solutions in 
this paper, to show global existence different methods will have to be applied. Function 
counting, i.e. the identification of free functions in the formal series solutions indicates, 
however, that spacetimes with a helical Killing vector and two regular horizons could 
exist: the series in the vicinity of the horizon contains two free functions which could 
be fixed in a way to allow for two regular horizons. Such an approach is suitable for a 
numerical treatment which could give - if successful - a strong indication of the existence 
for corresponding solutions to the Einstein equations. Though formal expansions do 
not provide existence proofs for solutions, they proved to be very useful and reliable to 
establish the behavior of cosmological solutions near singularities and served as a guide 
to prove existence and uniqueness, see |32i . 

To establish the Einstein equations in the presence of a helical Killing vector we use a 
projection formalism IH^ which leads to equations of 3-dimensional gravity with a 
sigma model as material source. The sigma model is determined by a generalized Ernst 
equation. The equations are discussed on the space of orbits of the Killing vector in a 
2 + 1-decomposition of the 3-space. The fixed points of the Killing vector, the horizons 
and the light cylinder where the Killing vector changes from being timelike to spacelike 
are singular points of the equations. They will be studied in terms of formal expansions of 
the metric functions in a chosen gauge. In this approach the horizons are regular horizons 
as in the case of the Kerr metric, but the series solutions contain two free functions which 
cannot be fixed locally. At the light cylinder the metric shows a similar behavior as at 
an ergosphere. If the constraints in the 2 + 1-approach are satisfied by the regularity 
conditions, the only equations to be solved are the two Ernst equations and the three 
'evolution' (with respect to the radial coordinate) equations. 

The behavior at infinity for spacetimes with a helical Killing vector is not yet well under- 
stood. There are indications from spherically symmetric models with equal amounts of 
ingoing and outgoing null dust jSH] ? which can be seen as a spherically symmetric analog 
to a helically symmetric spacetime, that the spacetime is not even asymptotically fiat 
in the weak sense that the metric tends asymptotically to the Minkowski metric. So far 
the only rigorous result due to Gibbons and Stewart [Hj is that there can be no smooth 
null infinity unless there is an additional axial Killing vector. Here we show that this 
is already case under the possibly too strong assumption that the spacetime tends to 
Minkowski spacetime asymptotically. Considering the Einstein equations in the presence 
of a helical Killing vector for large values of the radial coordinate, we show by using for- 
mal expansions of the metric that the spacetime is not asymptotically flat in the strong 
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sense that mass and angular momentum can be defined unless there is an additional axial 
Killing vector. The Weyl scalars behave asymptotically as 1/r, but the coefficients of the 
1/r terms have an oscillatory dependence on r. Thus there is no smooth X and no peeling 
in accordance with |6:. In asymptotically flat spacetimes for a stationary Killing vector, 
the Komar integral can be used to determine a conserved quantity via a surface integral 
calculated at finite radius, for which then the limit of an infinite radius is taken. We show 
that this limit is not defined under the above assumptions unless there is an additional 
axial Killing vector. 

The paper is organized as follows: In section 2 we use the projection formalism for the 
vacuum Einstein equations in the presence of a helical Killing vector. We discuss the 
resulting equations for the example of Minkowski spacetime. In section 3 we use 2 + 1- 
decompositions of the 3-space to study the singularities of the equations, the Killing 
horizons and the light cylinder. We use a formal expansion of the metric functions in 
local coordinates adapted to the singularities. In section 4 we study the linearized Einstein 
equations on a Minkowski background asymptotically in terms of formal expansions of 
the metric. We discuss the Weyl tensor and the Komar integral. In section 5 we add some 
concluding remarks. 

2 Quotient space metrics and Ernst equations 

The existence of a Killing vector can be used to establish a simplified version of the field 
equations by dividing out the group action. These quotient space metrics were ffist used 
in see also here we will follow |SZI. We use adapted coordinates in which the 
Killing vector ^ is given by ^ = dt where t is not necessarily a timelike coordinate. The 
norm of the Killing vector will be denoted by /. The decomposition we are using is 
not defined at the fixed points of the group action, i.e. the zeros of /, and the resulting 
equations will be singular at the set of zeros of /. The behavior of the solutions at these 
singular points will be discussed in the following section by studying formal expansions 
of the metric in the vicinity of the fixed points. 

In contrast to a standard 3 + 1-decomposition, the metric is written in this approach in 
the form 

I 

ds^ = -f{dt + kadx''){dt + hdx^) + -habdx^dx^] (1) 

latin indices always take the values 1, 2, 3 corresponding to the spatial coordinate. The 
Einstein equations in vacuum imply a Maxwell-type equation for what corresponds to the 
momentum constraint in a standard 3 + 1-decomposition (see |37j ) 

\Da{fV) = 0, (2) 

where kab = k^^a — ka,b {Da denotes the covariant derivative with respect to hab)- Notice 
that all indices here are raised and lowered with hab- If we define the twist potential b via 
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^^abc ^YiQ tensor density with e^'^^ = 1/ Vh) 

k"^ = -^—e^H „ (3) 

where h is the determinant of hab-, then equation is identically satisfied. The potentials 
/ and h can be combined to the complex Ernst potential £ = f + ib |HH1- The equations for 
/ and the integrability condition for b can then be combined to the generalized complex 
Ernst equation (the Ernst equation was originally obtained for the stationary axisymmet- 
ric case in [3^) 

fDaD'^S = DaSD'^S. (4) 

The 4 constraint equations in the standard 3 + 1-decomposition are thus replaced by a 
single scalar complex equation which is an advantage both for the analytical and the 
numerical treatment. 



The equations for the metric hab can be written in the form 

1 
2/ 



Rab = ^^{£,aB,b), (5) 



where Rab is the three-dimensional Ricci tensor corresponding to hab- Equations © 
describe three-dimensional gravitation with some matter model which turns out to be 
a 5*^(2, C)/S'0(1, 1) sigma model, see It is obvious that zeros of the norm of the 

Killing vector are singular points of the equations. 

To illustrate the above equations in the presence of a helical Killing vector, it is instructive 
to consider Minkowski spacetime in a rotating frame. In an asymptotically non-rotating 
frame, the Minkowski metric in standard cylindrical coordinates is in the above formalism 
for the stationary Killing vector given by / = 1, 6 = and hab = diag(l, 1, p^). In a 
rotating coordinate system where (p' = (p — Qt with constant Q, we get for a helical Killing 
vector = dt + Qdfp (notice that there is a helical Killing vector in Minkowski spacetime 
for arbitrary Q) 

/' = 1 - b' = -2nz, (6) 

where we have put a physically irrelevant constant in the definition of the twist potential 
equal to zero. Since the spatial metric hab in ((H) is rescaled by /, we have h'^f^ = f'hab 
except for h'^^ which is invariant under a transformation to a rotating frame. The light 
cylinder where the rotating observers corresponding to the vector ^ move with the velocity 
of light is given in this case hj p = In the interior of this cylinder, the Killing vector 
is timelike and / is thus positive, in the exterior it is spacelike. At the cylinder the 
signature of the metric hab changes from +3 to —1. In the four- dimensional picture there 
is no change in the signature of the metric but t and change roles at the light cylinder, 
(f) being a timelike coordinate in the exterior of the cylinder. 

The Ernst equation takes in non-rotating coordinates the simple form 

fAS = (Vf )^ (7) 
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where A and V are the standard differential operators in cyhndrical coordinates. In the 
rotating frame, the Laplace operator is replaced with the linear operator L defined by 

c- := dpp + -^dp + a,, + (1 - ^]V')^5<^<A, (8) 

which is just the helically reduced flat d'Alembert operator in a rotating frame. In the 
axisymmetric case (no 0-dependence) C reduces to the flat Laplace operator. In the non- 
axisymmetric case solutions to the equation C£ = behave for small p like solutions to the 
Laplace equation and for large p like solutions to a hyperbolic equation. Separating the 
angular dependence in spherical coordinates in a standard way via spherical harmonics, 
^ ~ J2ini^im{r)Yim{0,(l)), ouc recoguizcs that the solutions of CS = behave close to 
the origin as like solutions to the Laplace equation and for r ^ oo as e'^^^/r. One 
has thus to expect an oscillatory behavior for large r. Numerical studies of this type 
of equations have been carried out in [23 IHOl ISI] and [U E]- In general, Sommerfeld 
conditions (outgoing wave condition at finite values of r) have been used. Existence and 
uniqueness of solutions to boundary value problems for these equations were studied in 
using the theory of symmetric positive systems. 

In the stationary axisymmetric case the above equations can be further simplified. It can 
be shown (see e.g. [40j) that the spatial metric can be written in this case in the form 
hab = diag(e^'^, e^*^, p^). In this case the Ernst equation decouples from the equations for 
the metric function k and takes the form ((Tj) for a 0-independent Ernst potential. The 
metric function k follows for a given Ernst potential in terms of a line integral. In this 
formalism the Kerr solution for a single black hole with mass m and angular momentum 
J = im? sin ip takes a particularly simple form, 

e~*'^r+ + e*'^r_ — 2m cos 

O = : : , (9) 

e-*'^r+ + e*'^r_ + 2m cos 

where r± = ^J{z — m cos pY + P^i where the horizon is given by p = and \z\ < 
mcosif. 

Here we consider spacetimes with a single helical Killing vector. We adopt the definition 
of Friedman et al. Bl that the Killing vector can be written in the form ^ = dt' + Vtd^i 
where d^i is a spacelike vector with circular orbits of length 27r unless it vanishes. The 
vector dti is timelike outside the history of some sphere, is a constant, see jH] for details. 
This vector generalizes the helical Killing vector of Minkowski spacetime discussed above. 
It corresponds to the introduction of observers corotating with the binary system. Close 
to the black holes this vector will be timelike, but it will become null if these observers 
rotate with the velocity of light which determines the so-called light cylinder. 



3 2 + 1-decomposition, horizons and light cylinder 

In this work we are interested in spacetimes with a helical Killing vector that contain 
binary black holes. Since it was shown in ^ that spatially compact Killing horizons in 
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this case are event horizons, we are interested in spacetimes with two disconnected Kilhng 
horizons of spherical topology. 

With this assumption it seems convenient to introduce two systems of spherical coordi- 
nates adapted to the horizons in a way that one of them is given hj r = R = const. 
We assume that the spacetime can be globally foliated by spheres which is not necessary 
for the analysis below but for the planned numerical implementation. Since equations 
(jni) describe a model of three-dimensional gravity, it seems natural to use a 2 + 1 de- 
composition of the 3-space with respect to the radial coordinate. Let Afa = {A, 0, 0) and 
A/"" = (1, —B"')/A be the unit normal to the r = const surfaces; greek indices take the 
values 2 and 3. Denoting the metric of the r = const surfaces with Sa/3, we can write the 
metric hat in the form 

habdx^dx^ = + B''dr){dx'^ + B^dr) + A^dr^. (10) 

The Ricci tensor splits in the standard way (see e.g. 01]) in 3 'evolution' equations which 
contain second derivatives with respect to r, 

Rcp = ^iaMp) + 21C^^1C} - M^Mp + £mK^p + - K^pK (11) 

and 3 'constraint' equations below which have at most first derivatives with respect to r. 
Here ICajs is the exterior curvature of the r = const surfaces given by 

= ^V(a-B/3) - -S^p,-^ , (12) 

where Vq denotes the covariant derivative associated to Sap. The Lie derivative of Kap 
in the direction of Af is denoted by £f^fKa/3, 

1 ^7 /^7\ /^7\ 

and Ma = — (In^)^Q,. The two equations corresponding to the momentum constraint in 
the standard 3 -|- 1-decomposition read 

-RaaM"" = VpKi - Va/C^, (14) 

the Hamiltonian constraint is given by 

RabAf^AT^ - s^^Rap = 1C - ICaplC"^ - R^^^ . (15) 

where /C = /C^ and i?*^^-* = s"^R^^^. In the above equations, greek indices are raised and 
lowered with Sap and its inverse. 

For concrete calculations one has to fix a gauge. Since we are considering horizons of 
spherical topology and an adapted coordinate system, a natural choice of the metric Sap 
would be the standard metric of the 2-sphere. This is in general possible for smooth 
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metrics which we are interested in here since we are looking for regular horizons. Due to 
the fact that the spatial metric hat is rescaled by the norm of the Killing vector which 
vanishes at the horizon, the metric hab is expected to vanish there as well. A possible 
choice would then be a metric conformal to the metric of the 2-sphere, 

= /rMiag(l, sin^ 9). (16) 

This gauge will be convenient in the vicinity of the horizon. As the considerations for 
Minkowski spacetime in the previous section have shown, this choice is not possible across 
the light cylinder because of the signature change of the spatial metric there. Since we 
are interested in setting up equations suited for numerical treatment, we are looking for 
a system of coordinates which is able to cover the whole spacetime in the exterior of the 
horizon with a single coordinate patch. This could be possible with coordinates in which 
the angular part Qajj of the four-dimensional metric is the standard metric of the 2-sphere 
as in Schwarzschild coordinates. The latter would imply however the explicit inclusion of 
the vector ka in the equations which is contrary to the philosophy of the present approach 
to work only with its dual, the scalar twist potential b. 

Therefore we choose here a generalized Weyl gauge which we call 'quasi-isotropic'. We 
write 

s«^ = diag(rM^C). (17) 

A possible choice for C is C = r^sin^^^(l — R'^/r'^y. This corresponds to standard Weyl 
coordinates in which the horizon of a Kerr black hole is a sphere of radius R. These 
coordinates are related to the coordinates Q via p = {r — R^/r) sin 9 and z = {r + 
R? /r) cos9. The quasi-isotropic gauge thus reduces to the standard Weyl coordinates in 
the axisymmetric case. It can be used in principle throughout the light cylinder, but it 
remains to be shown whether it can be used globally. Due to the divergence structure of 
the Ernst equation which can be written free of covariant derivatives in the form 

f{h'''VhS,a),b = h'^'VhS^aSb, (18) 

the equation has in this gauge the standard terms of the Laplace operator for the rr and 
99 derivatives. This helps in the numerical treatment of the equations since standardized 
differential operators can be numerically inverted. In this gauge we have the non-vanishing 
Christoffel symbols corresponding to Sajj 

Tl, = i\nA),e,rl, = i\nA),^,Tl, = --1(1 - i?V^')' sin^ cos^, (19) 



and 



[1- Ryr^y sin' 9'' 



rL = „,.jt..2. ^r^3 = cot^. (20) 



(2) 

The components of the Ricci tensor read i?23 — ^^"^ 

= ■ = ■ - 1 - cot^(ln^),,. (21) 

(1 — it^/r^)"' sm fc' (1 — it^/r^)"' sm 
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The equations for Rat can be treated as in the case of a 3 + 1-decomposition: if the 
constraints are satisfied for some value of r, this will be the case for solutions to the 
evolution equations for all values of r. Since the horizon is a singularity for the equations, 
one has to give boundary conditions there which are compatible with the constraints 
and the evolution equations. With these boundary conditions one has to solve the Ernst 
equation which corresponds to two real equations and the three evolution equations (fTT|l ^. 
Thus one has to solve in total 5 equations as in the case of the IWM problem in |22j . The 
difference is here that the equations are not elliptic in the exterior of the light cylinder in 
contrast to the IWM equations and that the spacetime will not be asymptotically flat as 
discussed in the following section. 

To study the behavior of the metric at the horizon, we use formal power series in the local 
coordinate y = r — R. As in the case of ordinary differential equations, we adopt for a 
function F{r, 6, 0) the ansatz 

oo 

F(r,^,0) = l/"^^F,(^,0)y^ (22) 

j=0 

but here with coefficients Fj{6,(f)) depending on 6 and 0. The question is whether there 
are formal solutions to the Einstein equations of this form with vanishing norm / of the 
Killing vector for ?/ = which are more general than the Kerr solution for a single black 
hole. Here we are only interested in providing formal solutions intended for the use in the 
numerical treatment. Therefore we do neither discuss the convergence of the series nor 
global questions. We get (again we ignore a physically irrelevant constant in the definition 
of b) 



Proposition 3.1. In the gauge |7^, the equations ^ and (0j have formal solutions of 
the form / f^) with 



f = fo{e,<p)y'+fi{OA)y'+f2{0,<P)y'+---, b = bo{e,<p)y'+h{e,<p)y'+h{e,<p)y'+..., 

(23) 

and 

A = Ao{e, <P)y+Aiie, <P)y'+A2ie, (/))y^+. . . , = Qoie, <p)y^+. . . , 63 = $ o(^, <P)y'+. • • • 

(24) 

The functions /o(6',0) and bo{6,(j)) are free functions of 6 and with /o,(/,(O,0) = 0. All 
other coefficient functions in the expansions \2'Jji and \24^ can he expressed in dependence 
of fQ{6,(j)) and bQ{6,(j)), the leading order terms being 

A(^,0) = «:/o(^,0), ^ = -^' ^^ = -f' ^^ = "^' (2^) 

where the constant k, is given by k, = 2/(i?/o(0, 0)). 

^As in |42| . a fully constraint approach can be used alternatively in the sense that only the constraint 
equations instead of the evolution equations are solved. 



10 



The constant k in the relation between Aq and /o indicates a freedom in the choice of 
/o(0, 0). This freedom is due to the fact that a scale in the norm of the Killing vector is 
not fixed, after multiplication with some constant, ^ is still a Killing vector. 

Proof: 

With (jini) and we get for the exterior curvature ^ by using ^ and ^ 



/C22 = -Ao{R^ + 3Ry + 2y^)-Ai{2R^y + 5Ry^)-3R^A2y^ 

^y^^^Ao^ 
Ao 4 sin^ e 



+(60, - (In A),.eo)^ + ^^^y^ + O(y^) 



/C23 = Qeo,^ - (In A),,^eo^ + Q$o,e - cote^o^ ^ + 0{y^), 

■u^ 4 sin 6* COS 6*^ 7 / 'iin / n 

+*o,0;^ + ^,^3 Qoy^ + 0(2/^). (26) 

It is straight forward to check that the Hamiltonian constraint is satisfied to leading 
order (which is l/y^)- The momentum constraint ()14|) leads in lowest order to 

(lnA),„ = (ln/o),„. (27) 

Thus we have ^0 = ^/o with k = const. 

To ensure a regular axis {6 = 0) in spherical coordinates, the axis must be 'elementary 
flat', i.e. small circles around the axis must have an invariant circumference of 27r times 
the invariant radius in the limit of vanishing radius. This means for p = r sin 6, z = r cos 6 

lim / J'glldcj) = 27i lim 



\l9pp{p',z,(t))dp'. (2^ 
With the above relations this implies {ka is bounded at the horizon) 



Expanding this relation in y and using ()2H|1 . ()24|1 and ipTjl . we get in lowest order of y the 
condition that /o(0, 0) must be independent of and thus be a constant. This constant 
is related to n via 

The Ernst equation (jH) is satisfied in leading order for arbitrary /o(6',0) and bo{9,(j)), in 
the next higher order the real part implies 

/i + ^ = 0, (31) 

11 



whereas the imaginary part leads to 

-8&o/i + |&o/o + 56i/o = 0. (32) 

This determines /i and bi in dependence of /o and 60 • Iii order y^, the leading terms in 
the real part of the Ernst equation are 

(n - 2)2/o/n-2 (33) 

and 

n{n - A)hh^_^ + 8(2 - n)hofn-2 (34) 

for the imaginary part. The Ernst equation can thus be used to all orders determine /„_2 
and 6„_4 in dependence of quantities of lower order. 

With (j^mi we get for the Lie derivative of the exterior curvature (fT^ 



£j^lC^, = -tl3R + 2R'^]-A-7R^ + 2R'^-6R'^ 
y \ AoJ Ao A'q Aq 



40o,e . , ©0 , 3i?2(lnA),^$o 

--^-7(lnA).^ + ^-^ 



£nJC23 = ^(2eo,^-5(lnA),^eo + 2<l>o,e-4cot^$o-(ln^o),e$o) + 0(i/) 



-^-^^33 = ^F2— ( ^ + ^ ) JT- ( ^ + WWT- + ^ - ^ ) (35) 



9^1 SAj 2A2 
R^ ' 2RAo ^ T4f" A 

+ 16 sin cos 0^ -Asm'e{\nAo),e^, + ^ - 3(ln A),^^ + 0(i/) 



4 sin^ e 


(-- 








{2R 


AoJ 


Al [ 



Since there are second order derivatives with respect to r in the the evolution equations 
as in the Ernst equation, higher order terms in the expansion of the metric functions will 
appear here before they do in the constraints. Therefore there are no further conditions 
on the lowest order terms. In order 1/y the equation for R22 (there are no contributions 
from the Ernst potential in this order) leads to 

-(4 + 4i)=0. (36) 



y y^R Ao 

This implies that there are no terms of order y in the exterior curvature and no terms 
of order 1/y in the Lie-derivatives (jSSI)- 

The leading terms in order in the evolution equations are for n > 2 in R22 
-R\l - nf-^ + ^2 {nQn-3fi -{n + 3)(ln A),e©„_3 + cot ^e„_3) 

J^2 

(•l-^.s,,/,-!- (n- l)(lnA),^<l>„_3). (37) 



4 sin^ e 
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Similarly we get for R23 

{{n - l)en-3,ct> - 2n{lnAo),^Qn-3 + {n- l)$n-3,e - 2{n - 1) cot ^<l>„_3 - 2(ln A),e$„-3) 

(38) 

and for R^^ 

4 sin^ 6 

^2^4 (e„_3,e - 2(lnA),e0n-3 + ncot^e„_3) 
+ ^(n$„_3,^-2(ln^o),^'^'n-3)- (39) 

It is straight forward to solve the equations I38|l and for Qn-3 and $n-3- Function 
^„_i then follows from equation ^7\i . The Ernst equations (jH^ and (jMj) determine 
consequently /„_2 and 6„_4. We have thus shown that the evolution equations and the 
Ernst equation can be solved in this way to all orders. This completes the proof. 

The fact that An-i does not appear in the equations (jHSj) and implies that 6 and 
$ can be chosen to vanish for 0-independent /o and 60 • This leads as expected to the 
Kerr solution. Note that /„_2, &n-4 and An-i are determined algebraically by the above 
equations. Just to determine Qn-i and one has to integrate which leads to free 

integration functions. The latter are related to the fact that the used gauge conditions 
do not fix the gauge completely. There are transformations of the form 

r' = r + y^P{9,^) + 9' = 9 + y^S{9, 4>) + . . . , 4>' = 4> + y^T{9, 4>) + . . . (40) 

for non-trivial P, S and T which do not change the gauge. Since h'^.^ = hre + heg2yS and 
similarly for hr^, there are gauge modes in and $ which show up in the form of free 
integration functions. 

Remark 3.1. Due to the homogeneity of the Ernst equation in the Ernst potential, the 
functions fo{9, (p) and bo{9, 0) are not determined in the above expansions in the vicinity 
of a horizon. This gives hope that there might be a second regular horizon of spherical 
topology in the spacetime for a suitable choice of these functions. Whereas the behavior 
of the Ernst potential with respect to y is the same as in the case of a single Kerr black 
hole, the functions /o and 60 may be different. 

To treat the light cylinder we use a similar approach as for the horizon. By the definition of 
the Killing vector of |H] we are using here, the light cylinder will have cylindrical topology. 
We assume that the spacetime can be foliated by cylindrical surfaces, and use cylindrical 
coordinates in which the light cylinder is given by p = po = const. In an abuse of notation 
we use the same symbols for the 2 + 1-decomposition as used for the spherical case, 

habdx^dx^ = s«^(dx" + B^'dp) {dx^ + Bl^dp) + A^dr"^. (41) 

where dx"^ = dz. We use again the quasi-isotropic gauge which reads in this case 

Sa(3 = dmg{A\C). (42) 
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The choice C = is possible near the cyhnder. 

We assume that / has a zero of first order in v = p — po at the cyhnder since the Kilhng 
vector is supposed to change from timehke to spacehke there. Again we consider formal 
expansions of the metric function of the form 

oo 

F(p,z,0) = ^;'^-^F,(z,0y, (43) 

j=0 

We get 



Proposition 3.2. In the gauge the equations ^ and (0) have formal power series 
solutions of the form (|73| ) with 



f = f,{z)v + f,{z,(j)y + ..., b = bo{z) + b,{z,(j)y + ..., 

A = Ao{z)V^ + Ai{z,(P)v'^ + ... B2 = Zo{z,4>y + Bs = MzA)v- (44) 



The functions bo{z), Ao{z) and f2{z,(f)) are free functions of z and z, respectively. All 
other coefficient functions in the expansion ^4\ l can be expressed in dependence of bo, Aq 
and f2, the leading term being 

bt = fl (45) 

Proof: 

The formulas for the 2 + 1-decomposition in section 2 apply with the trivial change that 
r has to be replaced by p. The chosen gauge is also very similar to the one used for the 

horizon, the only difference being the factor in the expression for 822- Therefore we 

(2) 

will not give explicit formulas for the Christoffel symbols and i?)^^ here. For the exterior 
curvature (fT^. we get with 



JC22 = -^-UlV^+0{vl), 



/C23 = 7^$o,. + 0(t;i), 
/C33 = -T^fl-X^V^<^o,^ + 0(^^). (46) 

Aq^V \ Aq ) Aq 

This implies for the Lie derivative of the exterior curvature 



i^AfK23 = -^j-to,. + 0('t>°), 
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In lowest order, the equation for R22 then yields ()45|) . whereas the relations for R23 and 
i?33 can only be satisfied in this case for 

&o,^ = 0. (48) 
The 2;-component of the momentum constraint implies 

26o,.&2 + hzfo - ^hl, = 0, (49) 

which determines 62, whereas the 0-component of the momentum constraint gives 

= 0. (50) 

With condition P5|) the Hamiltonian constraint is satisfied to leading order. 

The real part of the Ernst equation gives no additional condition in leading order, in order 
V it leads to 

/i = (51) 

The imaginary part gives with in leading order (j^Uj) . Consequently 62 and /i are 
determined in this order. 

The equation for R22 reads with (fTT|l in order 1/v 

In the same order i?23 is identically satisfied. i?33 takes the form 

= V + -"^o,,/,. (53) 
>*-o Po 

This fixes $0 and Ai. Thus the functions &o(-2) and Aq{z) are not determined by the 
above equations. In addition Zq is not yet fixed. 

In higher orders of the expansion, the reasoning will be similar since the general structure 
of the equations is the same: the Ernst equation in order t>"~^ with n > 1 contains the 
leading terms n{n — 3)/o/„-i and 

n{n - 3)6„ + 26o,.%^ (54) 

in the real and the imaginary part respectively. The evolution equations in order t>"~^ 
contain the leading terms 



"'^-^-T -r + -72 Zn-2,z - {u - 2) (In A) ,^Z„_2 , (55) 
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and 

An 1 
(2n - 1)-^ + — - -j-2Zn-2,z- (57) 

Thus one can determine /^-i, &n, ^n, ^n-i and $n-i from the above equations unless 
n = 3. In this case the real part of the Ernst equation determines the function Zq which 
was still free, the function /2 remains undetermined. The equations in higher order fix 
all expansion functions except Ao, bo and f2- Free functions in Z and $ occurring after 
integration are again related to residual gauge freedoms as was the case in the vicinity 
of the horizon. Thus the series ()44j) provide a formal solution in the vicinity of the light 
cylinder. This completes the proof. 

Remark 3.2. It should be possible to apply Fuchsian methods as in j^SI to prove existence 
of the above solutions near the horizon and the light cylinder for some non- vanishing radius 
of convergence. However this would not answer the decisive question whether there can 
be two smooth horizons and a smooth light cylinder in the spacetime. Therefore we will 
not apply these methods here. 



4 Asymptotic behavior 

The formal solutions in terms of a series in the vicinity of the two horizons and the light 
cylinder in the previous section obviously do not imply global existence of a solution 
describing a spacetime with a helical Killing vector and two regular Killing horizons. The 
radius of convergence of these series is unknown. Therefore it is also not possible to 
make precise statements on the asymptotic behavior of the metric. As shown in [Hj such 
spacetimes cannot be asymptotically flat in the strong sense that they have a smooth 
null infinity. However this does not exclude the possibility that the spacetime is weakly 
asymptotically fiat in the sense that the spacetime tends to the Minkowski spacetime. This 
is what we will show in this section though the assumption of an asymptotic Minkowski 
metric might be too strong as indicated by the work of . However the used techniques 
will also be applicable for even weaker asymptotic conditions. 

Since the coordinate system we were using in the previous sections is asymptotically 
rotating (see the considerations for Minkowski spacetime in section 2), the Ernst potential 
is expected to have the kinematic terms (jH]). These terms will lead to technical difficulties 
if one wants to consider an expansion of the metric functions in powers of 1/r. Therefore 
we will consider in this section asymptotically non-rotating coordinates (t',r, 6',0') where 
0' = + Qt, t' = t. The metric will be studied via the linearized Einstein equations on a 
Minkowski background. Due to the helical symmetry the metric functions depend on t' 
and (f)' only via the combination x = (f)' — Qt'. 
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We assume the metric to be of the form qab = Vab + ^ab (capital indices take the 
values 0,1,2,3) for r — >■ oo where tjab = diag(— 1, 1, r^, sin^ ^) is the Minkowski metric 
in spherical coordinates, and where 6ab gives the deviation from Minkowski spacetime 
for large r. In cartesian coordinates 6ab is assumed to be of order 1/r. An analysis 
of the equations as in section 2 indicates however that the terms of order 1/r will have 
an oscillatory dependence in r (we exclude here possible logarithmic terms in the metric 
function). We will therefore consider a formal expansion of the metric functions of the 
form 

nrA4.)-t^5^, (58) 

j=0 

where the r-dependence of the Fj is to be understood to be purely oscillatory. In spherical 
coordinates we expect for the algebraic dependence on r that 6oo, 5oi and 6u are of order 
1/r, that 6o2, ^03, ^12 and 613 are of order r°, and that 522, ^23 and ^33 are of order r. 
We note that this ansatz for 5o3 allows for a so-called Newman-Unti-Tamburini parameter 
which corresponds to a magnetic monopole in electrodynamics. It is unclear whether such 
terms have to be expected in the present context, and the used methods are not suited 
to answer this question. 

To fix the gauge freedom we consider coordinate transformations of the form below which 
do not change 7]ab to leading order. The wanted gauge transformations can be put into 
the form 

r(r, 9, x) + ax 
t = t+ \, , 

, , , R{r,e,x) 

r = r + (j -\ , 

r 

Q{r,6,x) +'yx 



e + 



~ , , P{r,e,x) 

where the r- and x-dependence of the potentials T, R, Q and P is understood as before 
to be purely oscillatory, and where a, j3 and 7 are functions of 6 only. This implies for 
the transformed Minkowski metric in leading order 

^00 = -1 + — (Tx + a), goi = --{Tr + ^R,x), 5'02 = -^(Q,x + 7), 

_ 2 _ 

go3 = -nsin^ 9P^x, gii = l + -R,r, gi2 = Q,r, (60) 

r 

2/1 , \ ~ _ r-, , „;„2/io ~ .2„- 2/i/i , 2/9 



gi3 = sin 9P^r, g22 = r |^1 + —J , 5f23 = Q,x + sin 9P^e, ^33 =r sin 9 [1 

For metrics obtained by solving the linearized Einstein equations on a Minkowski back- 
ground, this implies that the functions T to P in (jHUI can be used to establish a certain 



17 



gauge. By an appropriate choice of P, we can choose ^23 to vanish under the made as- 
sumptions. By considering higher orders of this expansion, this should be possible to all 
orders. Similarly by choosing R and (3 we obtain ^22 = "^"^Qii- The additional freedom 
can be used to have a vanishing 5oo and (702- We assume that the gauge potentials T, R, 
Q and P are of the form T = EmGz(^m(^)e™^'^+^^^ + T-(^)e^'"(^-^'')). Therefore we had 
to add the functions a, [3 and 7 depending only on 6 to make sure that terms in the 5ab 
which are constant with respect to x can be compensated. This would lead, however, to 
terms proportional to x in 5qq and 5q2- Therefore we allow for a purely 6'-dependent 5qq 
and 5o2 which will not enter the linearized field equations to retain periodic potentials in 
X and r. Dropping the tilde, we thus choose the gauge (which is not exactly equivalent 
one used in the previous section, but close to it) 

^00 = -[1 + ^^^^ , 901 = ^, 903 = a, gss = sin^ 9 (^1 + , (61) 

and 

1 A Q 

9o2 = ho{e), gu = —g22 = 1 + — , 912 = — , 5-13 = ^- (62) 

The gauge is fixed up to a free function of r, in a, and functions of only in c, 9, and 
$, since we assume periodicity in r and x of the terms Fn{r, 6, x) in (j3Hj) . There could be 
a contribution in order to 0, but this must be a function of 6 alone as a consequence 
of the Einstein equations below and the periodicity condition in r. We put this function 
equal to zero here to fix a gauge freedom. This implies in leading order for the inverse 
metric 

g- = -(l-f^), = = <,- = ^(l-:^y (63) 

\ r J r r-^ sm 6 r-^ sm fc' V r J 



and 



02 ^0 11 2 22 1 ^ 13 ^ 12 ® /rA \ 

9 9 =rg =1 , g = — . 2 , g = — ^. 64 

r sm 6 



r 



We obtain 



Proposition 4.1. The linearized Einstein equations on a Minkowski background for the 
metric lead to two wave equations for the functions A and a, 

A^rr — fi^^, = 0, a^rr — fi^Ct,</)(^ = 0. (65) 

The remaining metric potentials follow in terms of quadratures. 
Proof: 

With relation ^ 

Rabcd = 7:{9ad,bc + 9bc,ad — 9ac,bd — 9bd,ac) (66) 
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for the linearized Riemann tensor we get for the Ricci tensor in lowest order in 1/r 



2 

r sin 



2ri?oo 


— — 2r2c ~ 2Vl? A ij)A) — VL^F A)A) 




2ri?oi 


= n{A + F),r<p 




2rRo2 






2-Ro3 






2ri?n 






2ri?i2 






2i?i3 






r 


— o'^Aj.j^ — A 




2-R23 






^2 ^^33 




(67) 



It is a consequence of the equations for i?03 and i?i3 that 

(]$,^ + a, = G'i(^), (68) 

where Gi is a free function of 6' only which is gauge invariant under transformations of 
the form (jSHI)- The equation for i?23 then implies 

^^a^^^ - a^rr = G2{r, 0), (69) 

where G2 is a free function of r and which reflects a gauge freedom and can be put 
equal to zero. Equation (j69|l represents the first of the two wave equations. We write the 
solution in the form of a Fourier series 

a = J2 e'""^{aiie)e'^'^'' + a" (^)e-'^'"'^). (70) 

meZ 

The reality condition for a implies ^. Thus we get for $ 

$ = - ^ 6*^^(0+ (^)e*^™" - a;^{e)e-'^''"') + Go{e). (71) 

The equation for Rqi implies 

A + F = G,{9,4>) + Ge{9,r). (72) 
It is then a consequence of -R22 (-R33 is identically satisfied) that 

A=J2 e""^A+ (^)e'^™^^ + A- (^)e-'^'"0. (73) 
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This gives the second wave equation. Again reahty of A imphes A_j^ = A^. 
Equations -Rqo and Ru lead to 

c=~ r A^^dr + G,ie). (74) 

These equations also determine that the right-hand side of ()72|) is only a function of 6 if 
the periodicity in r and (p is taken into account. Equations R02 and R12 then imply 

e = -^ r A,edr + G,{e). (75) 

This completes the proof. 

Remark 4.1. If the black holes have equal 'mass', i.e. equal combination of mass and 
angular momentum which can be defined via the Komar integral below, the spacetime 
has an additional discrete symmetry, it is invariant in a suitably defined coordinate system 
under the transformation —0. This implies for (I7n|l and (|7!?jl A^ = A" and = a~ . 
In this case no additional boundary conditions at infinity need to be given. A Sommerfeld 
condition which is typically considered at finite radius would only allow trivial solutions 
in this example if imposed at infinity. 

Remark 4.2. In case the functions a and $ have leading terms of order 1/r, i.e. if there is 
no NUT-parameter, the following equations for the Ricci tensor (jU7|) change 

2ri?i3 = Vta-r^ + fic 00 - A^r<f> + ^'^^m^ 

2ri?23 = l^a 50 + + e,^0 - A (76) 
These equations again imply wave equation (j^Hjl for a and 

^ = -k ^ (^,rd(P+l [ A^^dr + G7ie,r). (77) 



The ansatz (jHH|) already implies that the ADM mass and additional asymptotic multipoles 
cannot be defined due to the oscillatory behavior of the metric functions. We will show 
that it is also not possible to use the Komar integral asymptotically in a standard way to 
define a conserved quantity. This integral can be used to relate a locally calculated mass 
to the ADM mass for an asymptotically flat spacetime with a stationary Killing vector. 
The idea is to evaluate a surface integral at finite radius R and then to take the limit 
R ^ oo. Basically one uses that -kdC, is an exact differential which means that one can 
apply Gauss' theorem. We get for an integration over a sphere with t = const, r = const 



^[A,B]g''^g'''V^d9d^ = 0. (7J 

s 
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To calculate the integral near the horizon we need the inverse of the 4- dimensional metric 



= + A.A;^ g"' = -fk% g^' = fh^\ (79) 

where spatial indices are raised and lowered with hab- In the quasi- isotropic gauge we get 
with the results of section 3 for the surface integral that only the term 

/ V^ded(j)g''g'%,^^ (80) 
Js 

contributes. This leads with ^0 = i^fo to 

1 r f An 

- / sinOdedcj)^-^ = —. (81) 

2 Js A K 

For a single black hole the constant k = l/(2m). The Komar integral is of course only 
defined up to a scaling of the helical Killing vector, see the remarks in jS]. 

To check whether the surface integral can be defined for r ^ 00, we determine the 
integral for r = tq where vq ^ R and study whether the limit vq 00 exists. This could 
be possible if one uses the periodicity of the functions in in the 0-integration over a 
complete period at a finite value of r. As we will show below, this will not be the case 
because of the bilinear terms in the integrand. Note that the Killing vector reads in the 
used coordinates ^ = dt' + Vtd^i. It is readily seen that the integral can only exist if 
Gi = 0, since the corresponding terms in the integrand are of order r^. If we assume that 
this is the case, we get for the surface integral 



I I rsm9d9d(p{n'^^F^ + an{2 + Fr))- (82) 
^ Js 



The integral can only exist if these terms vanish after integration with respect to (p since 
the integrand diverges as r. Writing the integrands as a Fourier series as we have done in 
the proof of proposition 4.1, we get after integration with respect to (f) that the integrand 
is proportional to 

a+ + a- -^nJ2 m{aiA-e''''-^ - a" A+e'^^^'"'^). (83) 

This expression can only vanish if a^A^ + a^A^ = for m > 0. In the equal mass case, 
this is only possible if either or A^ vanish. If a and $ have leading contributions in 
order 1/r as in remark 14.21 the terms of order 0(r°)in the surface integral are due to (jHUj) 
and are of the form 

E ^i^^^^^"'"''^' - A^A-e-'^^n- (84) 

meZ 

The integral can only exist if the terms A^A^ vanish for all m 7^ 0. Since there cannot 
be purely 'outgoing' or 'ingoing' waves in the case of a helical Killing vector, this con- 
dition will lead to the axisymmetric case. Thus the surface integral cannot be defined 
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asymptotically in the presence of a helical Killing vector unless there is in addition an 
asymptotically axial Killing vector. In this case the integral just gives the expected value 
M + VlJ, the combination of mass and angular momentum corresponding to a helical 
Killing vector. 

Gibbons and Stewart [Hj showed that periodic boundary conditions are incompatible with 
a smooth null infinity. This is in accordance with the ansatz (j3Hj) as can be seen from the 
following consideration: We define the standard null-tetrad of Minkowski spacetime, 

r = -L(s„ + a,). „« = _L(_a„ + a^). . _i_ (a, + . ,35) 

We define the Weyl scalars as in [43J (we can use the above tetrad since we are only 
considering a linearization on a Minkowski background) 

C3 = 2Cabcdm''t'k'i'', 

C4 = -Cabcdk^^ik^m^ + tH^), 

C5 = 2Cabcdk''k'^^t''. (86) 



Determining the components of the Riemann tensor for the asymptotic metric of propo- 
sition 4.1, we get for the Weyl scalars in leading order 

Ci = — — {Arr + 2nAr^ + ^'^Afj,^), 

C2 = C3 = C4 = 0, 

C5 = -^{Arr-2nAr^ + n^A^^). (87) 

Thus the Petrov type is N. The Weyl scalars vanish for r 00, but this limit is in general 
not defined for rCi, i = 1, . . . ,5 because of the oscillatory behavior of the metric functions. 
Thus in accordance with jSj, there is no smooth X and no peeling in this case even if we 
assume that the metric tends asymptotically to the Minkowski metric. 



5 Outlook 

In the previous sections we have given a set of equations describing binary black hole 
spacetimes with a helical Killing vector. The equations have regular singularities at the 
Killing horizons and the light cylinder, and a non-regular singularity at infinity. This 
leads to a set of five equations which could be useful for a numerical implementation. The 
equations appear to be well suited for the multi-domain spectral method used in |22I, see 
also It is straight forward to include regular singularities in the spectral formalism 
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in the adapted coordinates we used for the analytical discussion, since the formal expan- 
sions we were discussing is very close to the philosophy of a spectral expansion. The main 
difficulty from a numerical point of view seems to be the oscillatory behavior at infinity. 
Typically a cut-off at some finite radius is used, but it is unclear which boundary condi- 
tions have to be used there. A possibility would be to match the solution at some large 
radius to an analytical solution of the linearized Einstein equations where it is not yet 
clear on which background the equations can be linearized (the asymptotic form of the 
solutions is still an open question). The main problem will be in any case the numerical 
resolution of the oscillatory metric close to infinity. 

From a mathematical point of view the most interesting question is whether there exist 
solutions with two regular Killing horizons in a vacuum spacetime with a helical Killing 
symmetry. In this paper we have only considered formal expansions of the metric in the 
vicinity of the singularities. The fact that the solution close to the horizons contains 
two free functions of the angular variables gives hope that such solutions might exist 
globally, but this needs to be proven. In case such solutions exist, it would be interesting 
to obtain the precise asymptotic behavior, whether the metric tends to the Minkowski 
metric asymptotically, and whether a NUT parameter is needed. Numerical results could 
give hints on how to answer these mathematical questions. 

The physical relevance of the studied model is clearly to obtain fully relativistic values for 
the ISCO and to get initial data for numerical calculations of the last phase of the binary 
system. In a real physical situation, the helical symmetry will be only an approximate 
symmetry. Therefore it would be interesting to study perturbations of a spacetime with 
an exact helical Killing vector studied here. There have been activities in this direction: 
in the Killing symmetry holds only in a finite region of space and time, the spacetime 
is asymptotically matched to a wave- zone. An approximate Killing vector was considered 
in 1^. 
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